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1.  Introduction 


The  objective  of  this  technical  report  is  to  provide  data  about  the  use  of  open  source  software  tools 
to  create  scientifically  correct  technical  animations  from  scientific  data.  The  packages  used  were 
Blender  3D,  blender3d.org,  and  the  Visualization  Toolkit  (VTK),  www.vtk.org.  Television  and 
movie  industry  animation  packages,  such  as  Maya1,  Houdini2,  Lightwave3  or  3D  Studio  Max1,  are 
very  powerful  and  very  costly.  These  programs  are  very  good  at  importing,  animating,  and  render¬ 
ing  the  millions  of  polygons  created  by  large  scientific  data  sets.  My  objective  was  to  see  if  we 
could  employ  an  open-source  package  with  the  same  or  nearly  the  same  results  as  using  a 
professional  package. 


2.  Discussion 


The  first  place  that  I  started  when  reviewing  Blender3D  was  the  number  of  polygons  supported; 
at  present,  version  2.42  supports  as  many  as  2  billion  faces  per  scene.  A  face  can  be  comprised 
of  quadrilateral  or  triangular  polygons.  On  average,  a  scientific  data  set  may  encompass  500,000 
to  1,000,000  cells  or  polygons;  larger  data  sets  contain  well  over  1,000,000  and  approach 
1,000,000,000  cells  or  polygons.  The  data  that  I  am  using  have  at  most  750,000  polygons; 
therefore,  Blender3D  will  have  no  problems  with  the  data. 

At  the  very  heart  of  Blender3D  is  the  Python4  programming  language.  Python  is  a  high-level 
programming  language  designed  by  Guido  van  Rossum  in  1990.  It  is  a  dynamically  typed, 
object-oriented  language,  which  makes  it  similar  to  Perl5,  Ruby,  Scheme,  Smalltalk6,  and  Tel. 
Python  runs  on  Windows7,  Linux8  or  UNIX9,  Palm  Handhelds  and  is  the  primary  programming 
language  for  Nokia  mobile  phones.  It  was  developed  as  an  open-source  project,  under  the 
Python  Software  Foundation.  Its  versatility  allows  for  many  different  uses  in  many  different 
areas.  National  Aeronautics  and  Space  Administration  uses  it  to  streamline  communication  with 
Oracle’s  databases  for  space  shuttle  missions,  the  New  York  Stock  Exchange  uses  it  to  create 


1  Maya  and  3D  Studio  Max  are  registered  trademarks  of  Autodesk,  Inc. 
2Houdini  is  a  registered  trademark  of  Side  Effects  Software,  Inc. 

■j 

Lightwave  is  a  registered  trademark  of  NewTek,  Inc. 

4Python  is  a  trademark  of  Python  Software  Foundation. 

5Perl  is  a  trademark  of  The  Perl  Foundation. 

SmallTalk  is  a  trademark  of  Active  Information  Corporation. 
7Windows  is  a  trademark  of  Microsoft  Corporation. 

o 

Linux  is  a  registered  trademark  of  Linus  Torvalds. 

9UNIX  is  a  registered  trademark  of  The  Open  Group. 
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graphical  user  interfaces  (GUIs)  for  the  real-time  equity  trading  system,  and  Industrial  Light  and 
Magic  use  Python  to  “glue”  together  the  thousands  of  computers  and  hundreds  of  software 
components  in  its  computer  graphics  pipeline. 

By  default,  Blender3D  does  not  provide  an  interface  for  color-per-vertex  shading.  We  needed  to 
find  a  way  to  map  scalar  values  from  the  data  to  the  vertices,  or  nodes,  of  the  polygonal  data.  I 
was  able  to  accomplish  this  through  the  use  of  VTK  and  a  new  Python  module,  VTKB lender .py, 
written  at  the  University  of  Alberta.  VTK  is  an  open  source,  freely  available  software  system  for 
three-dimensional  (3-D)  computer  graphics,  image  processing,  and  visualization.  VTK  is  a  true 
visualization  program,  meaning  that  it  does  not  just  allow  for  simple  geometry  visualization  but 
for  a  wide  variety  of  algorithms  that  include  scalar,  vector,  tensor,  texture,  and  volumetric 
methods.  It  also  allows  for  advanced  modeling  techniques  such  as  implicit  modeling,  polygon 
reduction,  mesh  smoothing,  cutting,  contouring,  and  Delaunay  triangulation. 

The  VTKBlender  Python  module  was  written  at  the  University  of  Alberta's  Academic  Information 
and  Communication  Technology  center.  The  module  converts  vtkPolyDataMapper’s  to  Blender3D 
polygonal  meshes,  with  color-per-vertex  information,  and  converts  Blender3D  meshes  to 
vtkPolyData. 

VTKBlender.PolyDataMapperT oBlender(pmapper,  me=None) 

If  the  function  is  run  with  only  one  argument,  this  function  takes  a  vtkPolyDataMapper 
pmapper  and  returns  a  new  blender  mesh  with  the  converted  polydata.  A  second  optional 
argument,  which  takes  a  pre-existing  blender  mesh,  may  be  provided,  causing  the  existing 
mesh  to  be  overwritten  with  the  polydata.  Please  note  that  the  new  mesh  is  not  added  to  the 
scene  and  can  be  added  afterwards  via 

sc  =  Blender . Scene . GetCurrent ( ) 
ob  =  Blender . Ob j ect . New (' Mesh ' ) 
ob . link (me) 
sc . link (ob) 

The  reason  why  the  function  takes  a  vtkPolyDataMapper  object  as  an  argument  (instead  of  a 
vtkPolyData  object)  is  because  the  vtkPolyDataMapper  can  also  contain  a  look-up  table  to 
color  the  data,  in  which  case,  the  blender  mesh  will  have  vertex  colors  set  accordingly. 

VTKBlender.BlenderToPolyData(me) 

This  function  takes  a  blender  mesh  and  returns  a  vtkPolyData  object  that  contains  the 
geometry  contained  in  the  mesh. 

(http://www.ualberta.ca/CNS/RESEARCHWis/VTKBlender/index.html) 


Blender  3D  is  an  open-source  software  package  used  for  3-D  modeling  and  animation  developed 
as  an  in-house  application  by  the  Dutch  animation  studio  NeoGeo  and  Not  a  Number  Technologies 
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(NaN).  Blender  is  available  for  use  on  many  operating  systems,  including  FreeBSD10,  IRIX11, 
GNU/Linux,  Windows,  Mac  OSX,  Solaris12,  SkyOS,  and  MorphOS.  Recent  additions  to  Blender 
support  for  2  billion  faces  per  scene,  user  interface  enhancements  and  support  for  third  party  ray¬ 
tracing  packages  have  brought  its  capabilities  close  to  those  of  “high  end”  3-D  software,  such  as 
XSI,  3D  Studio  Max,  and  Maya,  all  professional  programs  used  by  Hollywood  studios  such  as 
Pixar. 

The  user  interface  for  Blender3D  is  highly  customizable  and  configurable  by  the  user.  A  typical 
“scripting”  setup  will  have  an  OpenGL13  showing  the  data,  a  text  editing  window  showing  the 
Python  script,  and  application  windows  at  the  bottom  of  the  screen  comprised  of  rendering 
options,  directory  infonnation,  animation  options  and  hierarchical  data. 

The  Python  program  shown  in  figure  1  imports,  animates,  and  renders  the  time-variant  data  on  a 
local  file  system. 

All  polygonal  objects  in  Blender3D  are  considered  meshes.  The  VTKBlender  module, 
“VTKBlender.py,”  provides  the  interface  between  the  VTK  data  and  Blender3D  creating  the 
polygonal  mesh  and  applying  the  color-per-vertex  information  “on  the  fly”. 

The  data  represented  here  are  of  an  experimental  armature  for  an  electromagnetic  weapon 
system.  The  rendered  data  show  the  electromagnetic  forces  in  the  X  axial  direction.  In  the 
completed  animation,  I  was  able  to  move  the  camera  around  the  data  showing  multiple  views  of 
the  data.  The  areas  of  interest  (and  showing  where  the  forces  are  the  greatest)  are  in  red.  The 
completed  animation  presents  a  complete  picture  of  the  performance  of  the  candidate  annature, 
presenting  the  color-per-vertex  scientific  data  in  a  manner  that  is  easily  understood  by  the 
viewing  audience.  Figures  2,  3,  and  4  show  time-variant  electromagnetic  data  rendered  in 
Blender3D. 


10FreeBSD  is  a  registered  trademark  of  The  FreeBSD  Foundation. 
1  'iRIX  is  a  trademark  of  Silicon  Graphics,  Inc.  (SGI). 

Solaris  is  a  trademark  of  Sun  Microsystems,  Inc. 

IT 

OpenGL  is  a  registered  trademark  of  SGI. 
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Figure  1.  Blender3D  user  interface. 


Figure  2.  Blender3D  rendered  image,  time  step  1. 
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Figure  3.  Blender3D  rendered  image,  time  step  300. 


Figure  4.  Blender3D  rendered  image,  time  step  500  (the  last  time  step  in  data). 


3.  Conclusions 


In  the  past,  I  have  used  animation  packages  such  as  Maya,  Lightwave,  3D  Studio  Max,  and 
Houdini  to  present  and  render  scientific  data.  All  these  animation  packages  are  tremendous  at 
what  they  do,  but  there  is  a  cost  associated  with  each  package  and  in  some  cases,  it  may  be  in  the 
range  of  tens  of  thousands  of  dollars  each  year.  Employing  an  open-source  software  package 
such  as  Blender3D  to  render  photographic  quality  images  and  technical  animations  at  no  cost 
would  save  ARL  and  the  Department  of  Defense  tens  of  thousands  of  dollars  each  year  in  fees 
associated  with  software  licenses  and  maintenance. 
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The  incorporation  of  a  Python  interpreter  at  the  very  core  of  Blender3D  provides  an  application 
programming  interface  to  the  data,  thus  providing  a  means  to  read  and  manipulate  the  data  in  a 
highly  scientific  manner. 
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Appendix  A.  Blender  Python  Script  Importing  and  Animating  Time  Variant 
Data 


# !  BPY 


Name :  ' ANIMATE  XDMF ' 

Blender:  241 
Group:  'Animation' 

import  os, sys, string 
import  Blender 

from  Blender . Scene  import  Render 
try: 

import  VTKBlender  as  VB 

except : 

print  'Problem  importing  VTKBlender  module...' 

import  vtk  as  V 
#from  Xdmf  import  * 

from  XFrameWork . io  import  XdmfVtk  as  X 


class  Data: 

iteration  =  0 
datadir  =  None 
name  =  None 
ext  =  None 
me  =  None 

ActiveScalar  =  'FX' 
location  =  [[0,0,0]] 

def  _ init _ (self) : 

"""reads  XDMF  files  from  directory""" 
pass 

def  ComputeAverageDistance (self ) : 
x_vector  =  255 
dx  =  x_vector  /  18.0 

print  'Computing  average  distance...',  dx 
x  =  0 

for  i  in  range (  18  ) : 

loc  =  [  x  +  dx,  0,  0] 
self . location . append (  loc  ) 
x  =  x  +  dx 
print  ' X :  ' ,  x 

del  x_vector 
del  dx 

def  ImportTimesteps (self ) : 

"""Checks  to  see  if  we  have  multiple  timesteps,  if  so  import  them  into  Blender" 
print  ' %d  files'  %  self . iteration 
if  self . iteration  !=  0: 

self . Xdmf DataFiles . sort  ( ) 

for  file  in  self . Xdmf DataFiles : 

path  =  self. datadir  +  '/'  +  file 

self . Import_Xdmf (  self . Xdmf DataFiles . index (  file  ),  path  ) 

def  Import_Xdmf (self ,  cntr,  xdmf_file) : 

print  'XDMF  file:  ',  xdmf_file 
Blender . Window . WaitCursor ( 1 ) 
reader  =  X . vtkXdmf Reader ( ) 
reader . SetFileName (  xdmf_file  ) 
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reader .Update Informat ion ( ) 
reader . EnableAllGrids ( ) 
reader . EnableAllArrays () 
reader . Update ( ) 

surface  =  V. vtkDataSetSurfaceFilter ( ) 
surface . Setlnput (  reader . GetOutput ( )  ) 
surface . SetPiecelnvariant (  1  ) 
surface . Update ( ) 

cells  =  surface . GetOutput (). GetNumberOf Cells ( ) 

lut  =  V. vtkLookupTable ( ) 
lut . SetNumberOf Colors (  256  ) 
lut . SetHueRange (  0.667,  0.0  ) 
lut.  Build  () 

m  =  V. vtkPolyDataMapper ( ) 
m. Setlnput  (  surf ace . GetOutput  ( )  ) 
m. SetScalarModeToUsePointData ( ) 

m. Getlnput ( ) . GetPointData ( ) . SetActiveScalars (  self . ActiveScalar  ) 
min,  max  =  m. Getlnput (). GetScalarRange ( ) 
if  min  <  0.0: 

min  =  0.0 

print  'Scalar  min:  %f  max:  %f  '  %  (  min,  max  ) 
m. SetLookupTable (  lut  ) 

m. GetLookupTable ( ) . SetTableRange (  min,  max  ) 
m. Update ( ) 

self . name, self . ext  =  os .path . splitext (  os .path.basename (  xdmf_file  )  ) 
sen  =  Blender . Scene . GetCurrent ( ) 

if  entr  ==  0: 

self. me  =  VB . PolyDataMapperToBlender (  m  ) 
ob  =  Blender . Ob j ect . New (' Mesh ' ,  self. name) 
ob . link (self .me)  #tell  object  to  use  new  mesh 
sen. link (ob)  #link  object  to  current  scene 

else : 

self.u_me  =  VB . PolyDataMapperToBlender (  m,  self. me  ) 
c_ob  =  scn.getActiveObject () 
sen. unlink  (  c_ob  ) 

ob  =  Blender . Ob j ect . New (' Mesh ' ,  self. name) 

ob . link (self ,u_me)  #tell  object  to  use  new  mesh 

sen. link (ob)  #link  object  to  current  scene 

cam  =  Blender .Object . Get (' Camera ' ) 

grpl  =  Blender .Object . Get (' armature_guides ' ) 

grp2  =  Blender .Object. Get ( 'payload' ) 

lampl  =  Blender .Object . Get (' Lamp ' ) 

lamp2  =  Blender . Ob j ect . Get (' Lamp . 001 ' ) 

ob .makeParent (  [cam,  grpl,  grp2,  lampl,  lamp2]  ) 

ob . setLocation (  self . location [entr]  ) 

ob.sel  =  1 

ob. Layers  =  sen. Layers 
sen .update ( ) 

Blender . Window . RedrawAll ( ) 
del  sen 

def  RenderScene (self ,  entr): 

context  =  Blender . Scene . GetCurrent (). getRenderingContext ( ) 
context.saveRenderedImageC%s.%04d.jpg"  %  (  self .ActiveScalar ,  entr)  ) 
context . setlmageType (Render . JPEG) 
context . enableExtensions ( 1 ) 
context . render ( ) 

del  context 

def  Main (self) : 

self . ComputeAverageDistance ( ) 
for  i  in  range (1) : 
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file  = 

'  /home/ j vines/Devel/Coupling_05/EM_DYNA3D_COUPLING/TEST/Coupled_PROBED_%04d . xmf ' 

self . Import_Xdmf (  i,  file  ) 
self . RenderScene (  i  ) 


if  _ name _ ==" _ main 

d  =  Data ( ) 
d.Main ( ) 
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